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Abstract. Generic non-magnetic disorder effects onto those topological quantum 
critical points (TQCP), which intervene the three-dimensional topological insulator 
and an ordinary insulator, are investigated. We first show that, in such 3-d TQCP, 
any backward scattering process mediated by the chemical-potential-type impurity is 
always canceled by its time-reversal (T-reversal) counter-process, because of the non- 
trivial Berry phase supported by these two processes in the momentum space. However, 
this cancellation can be generalized into only those backward scattering processes which 
conserve a certain internal degree of freedom, i.e. the parity density, while the 'absolute' 
stability of the TQCP against any non-magnetic disorders is required by the bulk-edge 
correspondence. Motivated by this, we further derive the self-consistent-Born phase 
diagram in the presence of generic non-magnetic disorder potentials and argue the 
behaviour of the quantum conductivity correction in such cases. The distinction and 
similarity between the case with only the chemical-potential-type disorder and that 
with the generic non-magnetic disorders are finally summarized. 
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1. Why non- magnetic disorders in 3-d quantum spin Hall systems ? 

A 3-d Z 2 quantum spin Hall insulator [U EJ |3] is denned to be accompanied by the integer 
number of surface conducting channels, each of which is described by the 2+1 massless 
Dirac fermion, i.e. helical surface state. The stability of each helical surface state is 
protected by the Kramers degeneracy at the time-reversal (T-) symmetric (surface) 
crystal momentum. The associated spin is directed within the XY plane, and rotates 
odd-number of times when the surface crystal momentum rotates once around this 
T-symmetric /c-point. While an individual helical surface state is protected by the T- 
symmetry, level-crossings between any two different helical surface states are generally 
accidental: they can be lifted by a certain T- symmetric perturbation [3l HJ EE E] . As 
a result of this pair annihilation, those insulators having even number of helical surface 
states can be adiabatically connected to an ordinary band insulator having no surface 
conducting channels at all ('weak topological insulator'). On the other hand, those with 
odd numbers of helical surface states are always accompanied by (at least) one gapless 
surface state, even when subjected under this pair annihilation ('strong topological 
insulator'). In the presence of the T-symmetry, the latter type of insulators cannot be 
decomposed into any two copies of spinless wavef unctions, and therefore regarded as a 
new quantum state of matter which goes beyond the quantum Hall paradigm. 

The quantum critical point [3 El El \10\ which intervenes this 3-d topological 
insulator and a T-symmetric ordinary band insulator is also non-trivial. That is, the 
stability of its critical nature is tightly connected to the stability of each helical surface 
state in the topological phase. To be specific, consider that a certain T-symmetric 
model-parameter is changed, such that the system transits from the topological phase 
to a T-symmetric ordinary phase. During this, the helical surface state in the topological 
phase is always protected by the T-symmetry. When the system reaches the quantum 
critical point, however, the bulk-wavefunction becomes extended once. Thus, the two 
helical surface states localized at two opposite sample boundaries can communicate via 
this extended bulk-wavefunction, only to annihilate with each other. Thanks to this 
pair annihilation mediated by the extended bulk, the system can safely remove the stable 
helical surface state and enter the ordinary phase having no surface conducting channels, 
while simultaneously keeping the T-symmetry. To put this reversely, the bulk-extended 
character of the quantum critical point is required to be stable against any T-symmetric 
perturbations, provided that an individual helical surface state in the topological phase 
is stable against the same perturbations. Such a quantum critical point can be dubbed 
as the strong topological quantum critical point. 

The above picture in the clean limit implies non-trivial disorder effects around 
this topological quantum critical point [UJ. To describe this with clarity, let us first 
define three independently-tunable physical parameters in the topological phase: (i) 
the (topological) band gap W topo , (ii) the band width of the conduction/valence band 
Hband, and (iii) non-magnetic disorder strengths W<n s - Provided that Wdis < W top o, the 
helical surface state in the topological phase is stable against any non-magnetic disorders 
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irrespective of Wband [121 ESI EH]- On increasing W^is such that WWd < W<n s < W topo , 
those bulk-states far from the two band-centers, i.e. that of the conduction band and of 
the valence band, become localized (see Fig. 1(a)) \IE\- Especially near the zero-energy 
region, /i = 0, the system has no extended bulk wavef unctions. Thus, two helical surface 
states localized at the two opposite boundaries respectively are disconnected from each 
other. Starting from such a localized phase, consider decreasing the topological band 
gap (or increasing the disorder strength), such that Wband < Wtopo < Wdis- During this 
decrease, these two disconnected surface states should communicate with each other 
once, so as to annihilate with each other, before the system enters the ordinary insulator 
phase. This requires that, at Wt op o — Wdis, the bulk-wavefunctions at \i ~ must 
become extended once, so as to mediate these two opposite boundaries. Moreover, it 
also requires that such a bulk-extended region always intervenes between the topological 
phase and the ordinary insulator phase in the three-dimensional parameter space 
spanned by the chemical potential fi, W<n s and W^opo (see Fig. l(b,c)). Such behaviour 
of the bulk extended region is often dubbed as the levitation and pair-annihilation 
phenomena, and was ubiquitously observed also in other topological systems, such as 
quantum Hall systems [161 LTZl HE] and 2-d Z 2 quantum spin Hall systems [191 EDI EI]- 




Figure 1. (a) A schematic density of state at W^is < W / t0 po(= rri). (b,c) schematic 
phase diagrams in the n-Wt opo -Wdi s plane. The blue regions in these three figures stand 
for the delocalized state, which always intervene between the topological insulator 
phase and an ordinary insulator phase in any 'T-symmetric' phase diagram. The 
figures are derived from Ref. [IT] . 

Now that the 3-d Z 2 topological band insulator is a new quantum state of matter 
having no quantum-Hall analogue, the bulk-delocalized phase which intervenes this 
insulating phase (or localized state adiabatically connected to this) and the ordinary 
band insulator (or corresponding localized state) is also a new type of three-dimensional 
metallic phase, which has no analogue in any other condensed matter paradigms. Taking 
into account the chemical-potential type disorder (the most representative non-magnetic 
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disorder), we have previously studied [TT] the self-consistent Born phase digram around 
this topological metallic phase (or quantum critical point) and calculated the weak- 
localization correction to the charge conductivity. In actual solids [2H ES, EEJ [27J EE] , 
however, electronic disorders can generally exist not just in the chemical potential, but 
also in the transfer integral and spin-orbit interaction potential itself. On the one hand, 
the above argument based on the bulk-edge correspondence clearly dictates that the 
critical nature of the topological quantum critical point is stable against any of these 
non- magnetic disorders (not just against the chemical potential type disorder). 

Motivated by this, we will treat in this paper various types of non-magnetic disorder 
potentials on an equal footing, so as to argue actual/generic properties of the topological 
quantum critical point in disordered Z 2 quantum spin Hall systems. The organization of 
this paper is as follows. In sec. 2, we summarize the effective continuum model studied 
in this paper. We next argue in sec. 3 that the non-trivial Berry phase (i.e. e m ) inherent 
in the 3-d topological quantum critical point (TQCP) induces the perfect cancellation in 
the backward scattering channels mediated by the chemical-potential type disorder. This 
partially upholds the above argument based on the bulk-edge correspondence. However, 
it also becomes clear that this 'Berry phase' argument does not necessarily work in the 
presence of those non-magnetic disorders other than the chemical-potential type one. 
More accurately, those backward scattering processes which do not conserve a certain 
internal degree of freedom, i.e. parity density degree of freedom, are not generally set 
off by their T-reversal counter processes. Promoted by this theoretical observation, we 
further derive in sec. 4 the self-consistent Born phase diagram in the presence of generic 
T-symmetric disorder potentials. Based on this, we argue in sec. 5 the weak-localization 
correction to the charge conductivity in the presence of both the chemical-potential type 
disorder and the topological-mass type disorder. The distinctions between the case with 
only the chemical potential disorder and that with generic T-symmetric disorders are 
summarized/clarified in sec. 6. 

2. effective continuum model and T-symmetric disorders 



An effective continuum model describing the topological quantum critical point is given 
by a following 3 + 1 Dirac fermion [31 El |9] ; 



where the two massive phases, m > and m < 0, represent the quantum spin 
Hall insulator and the T-symmetric conventional (band) insulator respectively: 'm' 
corresponds to the topological band gap Wtopo- The 4 by 4 7-matrices consist of the 
Pauli spin matrix part s and the other Pauli matrix a representing the sublattice (or 
orbital) degree of freedom. There exist at most five such 7-matrices which anti-commute 
with one another. Here they are dubbed as 71, 72, 73, 74 and 75. The product of these five 
always reduces to (the minus of) the unit matrix; 7172737475 = — 1. Thus, even numbers 




(1) 
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of 7-matrices out of these five should be time-reversal (T)-odd and spatial inversion (X)- 
odd. In eq. (pQ), three of them are already linearly coupled with the momentum. Thus, 
four of them should be both T-odd and X-odd. Without loss of generality, eq. ([1]) took 
7i, 2,3,4 to be T-odd and X-odd, so that 75 is T-even and X-even. Under this convention, 
7i5 = -*TiT5> 725 , 735 and 745 are all T-even and X-odd, while 712, 713, 7i4, 723, 734 
and 742 are all T-odd and X-even. Thus, arbitrary 'on-site' type non-magnetic disorder 
potentials are generally given by, 

r 4 

Himp = d 3 ri)\r){v Q (r)% + v 5 (r)% + ^v j5 (r)j j5 }'4)(r). (2) 
J 3=1 
with the real- valued functions Vj (j = 0, 5, 15, • • • , 45). 

In Ref. [TT], we have further assumed that the chemical-potential type disorders is 
the most dominant and have set the other five to be zero: v$ = V15 = ■ ■ ■ = t> 45 = 0. 
Based on these simplifications, we have carried out the self-consistent Born analysis and 
weak-localization calculations, and obtained a simple microscopic picture, which can 
indeed explain how the bulk-extended region emerges when the topological band gap 
'to' changes its sign [TT]. In the next section, we will describe an alternative argument, 
which more directly dictates that the topological metallic phase (or quantum critical 
point) is free from any localization effect induced by generic non-magnetic disorders. 



3. 'Berry phase' effect in 3-d topological metallic phases 

Our argument here is a straightforward 3-dimensional generalization of the 'Berry phase' 
argument invented by Ando et al. [22] in the context of the carbon-nanotube (or 2-d 
graphene) subjected under long-range impurity potentials. Specifically, we will prove 
that an individual backward scattering process is precisely canceled by its time-reversal 
counterpart, which leads to a complete absence of the backward scattering in the 3-d 
topological quantum critical point. 

Consider first the most simplest case, where the T-matrix is composed only by the 
chemical-potential- type impurity vo(r), 

* 1 11 

In the momentum representation, the (p + l)-th order backward scattering term reads, 

v (-k -hp)--- v (k 2 - h)v (ki - k) 



(-k,e- k ,a\f^\k,e k ,a')= £ ••• ]T 



(-k,e- k ,o-\k p ,T p e kp ,a p ) x ••• x (k 2 ,T 2 e k2 ,a 2 \k 1 ,T 1 e kl ,a 1 )(k 1 ,T 1 e kl ,a 1 \k,e k ,a'), (3) 

where Tj = ± (j = 1, ■ • • ,p) in combination with e kj specifies the eigen-energy of the 
hamiltonian at the momentum kj, i.e. Tje k . = ±|%|- Such an eigen-energy is always 
accompanied by doubly degenerate eigen-states, which is specified by the spin index, 
(Tj =t, |. Namely, each eigen-state at the momentum kj is uniquely identified by these 
two indices, \kj,Tje k .,aj). Without loss of generality, we can take the initial ket-state 
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and final bra-state to have an positive eigen-energy, i.e. \k,e k ,a) and (— k, e_fc, a\. We 
can also fix the gauge of the eigen-wavefunctions, e.g. 

\k,re k .,a) = e <§(^«*-«*^)|r,<7) 3 ee U k \r,a) 3 , 

with k = \k\(sm9 cos 0, sin^ sin 0, cos#). |r, (7)3 is the two-fold degenerate eigenstates 
of 73 which belong to its eigenvalue r, i.e. 73 |r, <t) 3 = r|r, cr} 3 . 

Based on this set-up, we can explicitly argue that the following two backward 
scattering processes, which are related by time-reversal operation, set off each other in 
eq. ©, 

\k, +, a) -> \h, n, ■ ■ •) ->■ ► r p , •••)-> | — fc, +, cr'), 

I A;, +, o") -+ I - k p ,T p , ■ ■ -} -»- ► I - fci, n, •••)-»• I - A;, +, a'). 

where the summation over the spin indices associated with the intermediate states are 
assumed. To be more explicit, we can show that the following two have an opposite sign 
with each other, 

S { ^ 1) = y^(-fc, e k , a\k p , r p e kp , a p ) ■ ■ ■ (h,Tie kl ,ai\k,e k ,a'), 
to} 

S iJ l) ' = ^2(~ k ^k,o-\ - &4,Tie fcl ,cxi) ■ • • (-k p ,T p e kp , a p \k, e k , a'). 

for any a and a'. To see this, take the sum over all the intermediate spin indices first. 
This leads to 

S^ 1} = 3 (+, o-\U[ k ■ (70 - np^Jn) ■■■(%- n 2 , u %) ■ (70 - ™i,a7a) -U k \+, o~') 3 . 

S^ 1} ' = 3 (+, a\U ] _ k ■ (70 + tii jM 7 m ) • (7o + n 2jV %) ■ ■ ■ (% + n p ^x) -U h \+, cr'} 3 . 

where rij is a product between the normalized vector parallel to kj and Tj, rij = Tjkj/\kj\. 
The sequential product of the projection operators can be further expanded in this unit 
vector, i.e. 

(1 + n p ■ 7) • • • (1 + n 2 • 7) • (! + n i ■ 7) 

= 1 + - 7 + ^ (nj •7)(n i -7) + • •• + (n p • 7) • • • (n 2 • 7) (m -7)- (4) 

*>i 

An even-order term in this expansion is always a linear combination of 70, 723, 731 and 
712, 

(n 2l ■ 7) • - • (ni • 7) = ^70 + e^xBi^yx. 

Ai is a scalar quantity, which is composed only by inner products among 21 unit vectors. 
On the other hand, Bi clearly behaves as a vector, being coupled to 723, 731 and 712. It 
is expressed as (a linear combination of) the product between I — 1 inner products and 
one outer product, 

Bi^= b ljM-no( n i xn j)^i n i-n m )---{n n -n ), 

l,m,---,n,oj^i,j 
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with some coefficients b'Z free from • • • , n 2 i}. Multiply this even-order term by one 
additional n ■ 7, one can obtain any odd-order term in eq. 01]) as a linear combination 
of 7b 72, 73 and 745; 

{n 2 i+i ■ 7) ■ • • (ni • 7) = Q^jf, + A745- 
C\ is a vector composed of I inner products, while Di is a scalar containing one scalar 
triple product, 

^ ^ C m,no - pq n rn,Lt( n n ' n o) ' ' ' ( n p ' n q)i 
n,o,---,p,q^m 

D i= Yl dl i jm ,no- Pq ( n i xn j)- n Unn-n )---( y n p -n q ). 

n,0—,p,q^i,j,m 

Thus, under the time-reversal operation, i.e. (n p , n p _i, ■ ■ ■ , ni) — > (— ni, —n 2 , ■ ■ ■ , — n p ), 
A\ and D\ do not change the sign, while B\ and C\ change their signs. These observations 
lead to a perfect cancellation between S^ 1 ^ and S^ 1 ^' , 

S^ 1] + = 2^2 3 (+, o-Wh ■ {A% + £745} • U k \+, oOa, 

= 2i^2 3(+,o-|(- sin 0723 + cos07 3 i) • (A% + Dj^ 5 )\+,a') 3 

a" 

for arbitrary cr and a'. This dictates that the backward-scattering process mediated by 
the chemical-potential-type disorder vanishes completely at the 3-d TQCP. 

More generally, we can prove that any backward scattering process which conserves 
the eigenvalue 0/745 is always set off by its corresponding T-reversal counterpart. Such 
a process is always mediated by even number of either Vjs(r) (j — 1, 2, 3) or v$(r), since 
both 7j 5 (j = 1,2,3) and % anticommute with 745. To uphold the complete absence 
of the backward scattering in this case, we have only to show that the following two 
matrix elements have an opposite sign with each other, 

S^ 1} = 3 (+, (y\Ul k la p {lo - rip^f,) ■ ■ ■ 7 a2 (7o - n 2 ,u%)% 1 (% ~ n 1>x A fx)% Uk\+, °~ ')s 

S a£ = 3<+, O"l^lfc7a (70 + n hjJi %)j ai (% + n 2 ,v%)% 2 • • - (70 + ^p,A7A)7a p ^fc| + , 0-/3, 

for arbitrary a and a'. The indices ao, ai, a 2 , • ■ ■ and a p can be either 0, 15, 25, 35 , 45 or 
5, under the condition that the total number of those 715, 725, 735 and % which mediate 
the initial state and the final state is always even. To see the relative sign between 
these two, let us first transport all the intervening 7 aj s in the leftward/rightward until 
they meet with the bra/ket states in S^ 1 ^ / S^ 1 ^' . This transport brings about an 
appropriate redefinition of the normalized vectors, rij — > Tij, 

s *£ = (y\U^ k la p ■ • -7ai7a • (to - ™p,m7m) " ' (% ~ ni,\l\)Ul\+, a') 3 

S^ 1 '' = 3 (+, (y\Ul k {lO + ni jjU 7 M ) • • • (70 + ™p,a7a) • laola, ■ ■■la p U k \+, J)*. 

S( p+1 ^ and S( p+1 ^' is still connected by the T-reversal operation, (n p , n p _i, ■ • • , Tii) <->■ 
(—Til, • • • j —Tip-i, —Ti p ). Thus, the preceding expansion can be used again, 

(70 - ^ p ,m7m) • • • (7o - ^Tv) (70 - ^i,a7a) = -47o + B^t^xlux + + D% 5 , 
(7o + ni )At 7 M )(7o + n 2)U %) • • • (70 + ^p,a7a) = -47o - B^xluX ~ C^Tu + D745. 
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The condition imposed on a,j dictates that 7 Qp • • •7 ai 7 ao and 7 ao 7ai • • 4 lap always 
reduce to either (i) and — 7, (j = 1, 2, 3, 12, 23, 31) respectively or (ii) ^ m and 7 m 
(m = 0,45) respectively. Consider the former case first, e.g. 

S£t 1} = 3 <+, <rp_ k ■ 7i • (Alo + B% x + C% + D% 5 ) -U k \+, a') 3 

= s<+, o-l^lfc^fc • m p7 p • (A 7o + 5'i A + C 7 '^ + D% 5 )\+, a'} 3 , 

S^' = o-jC/ljfe • (-A70 + B% x + C 7m - ^745) • 7i • &k\+, cr') 3 

= 3 (+, a\JJ ] _Jj k ■ (-A70 + £7„a + C 7 P - D745) ■ m p 7 p |+, o-'} 3 . 

with U ] _Jj k = - sin 0723 + cos07 3 i and m p j p = U^iUk- When % p ■ ■■% 2 % 1 = 712 (or 
723, 731), r«p7p should be replaced by m p e p ^7^. In either cases, S^t and S^t set 
off each other completely e.g. 

S^ 1] + S^ l) ' = 3 (+, a|(-^ 7 23 + W3i) • {5' 7 ,a + ^7m. ™p7p>I+. Ob 

= 2 x 3 (+,o r |(-^723 + C073i) • (^m^ + (^^70)1+, cr') 3 
= (---)x 3 (-,a / |+,a / ) 3 = 0. 

A similar cancellation also holds true for the case with (ii) ^ ap ■ ■ ■ ^ a2 la 1 = 1^1 a 2 4 4 4 la p = 
745. These observations thus conclude that any backward scattering process which 
conserves the eigenvalue of 745 is always canceled by its time-reversal counter-process. 
Like in the 2-d graphene case [22], this cancellation is actually a direct consequence 
of the non-trivial Berry phase 'e i7r ' accumulated along a closed loop composed by two 
time-reversal counterparts in the momentum space (see Fig. [2]). 




Figure 2. Red curves denote two backward scattering processes on a sphere in the 
momentum space (\k\\ = ■ ■ ■ = \k p \ = \k\), which are mutually time-reversal to each 
other. The blue-shaded region, whose boundary is shaped by this two path, occupies 
one half of the whole surface of a sphere. Thus, an electron acquires the Berry phase 
7r (instead of 2tt), when it travels once around this boundary. 

The above 'Berry phase' argument, however, does not hold true for the case with 
those backward scattering processes which change the eigenvalues of 745. To see this, 
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consider, for example, the backward scattering processes mediated by odd number of 
f 5 (r) and compare following two T-reversal paired processes, 



We have already made any two 'neighboring' 75 to annihilate each other, which leads 
an appropriate redefinition of the normalized unit vectors. The extra 75 was then 
transported leftward/rightward in S^ 1 ' / S^ 1 '' , until it meets with the bra/ket-state. 
Exploiting the expansion described above, we can readily see that these two do not 
cancel exactly in general, 



Stt l) + S^ l) ' = 2iJ2 3<+, <r\(- sin + cos0 7 3i) ■ (A + C 3 %) ■ %\+, a% ? 0. 



under short-range impurity potentials, where the intervalley scatteringjj suppresses the 
cancellation due to the Berry phase and consequently induces the crossover from the 
symplectic class to the orthogonal class j22j EH]. The situation in the 3-d TQCP is, 
however, quite different from this 2-d analogue. First of all, the preceding argument 
based on the bulk-edge correspondence requires the extended character of the 3-d TQCP 
to be stable against any non-magnetic disorders. Moreover, it always belongs to the 
symplectic class irrespectively of the type of non-magnetic disorders, which was not 
the case in the 2-d graphene [23]. Therefore, not only from the actual material's 
viewpoint [2U [25j [26j |27J [28], but also from the theoretical standpoint, it is quite 
intriguing and important to investigate the effect of general T-symmetric disorder 
potentials (including the topological-mass type one) in the 3-d topological metallic 
phase. In the next section, we will first study how a single-body Green function is 
renormalized by these general non-magnetic disorders. Based on this analysis, we will 
argue the behaviour of the quantum conductivity correction in the presence of both the 
chemical-potential type disorder t> and topological-mass type disorder w 5 . 

4. self-consistent Born phase diagram — general case — 

Consider Hq + "Hi mp with randomly distributed real-valued Vj(r) (j = 0,5,---,45). 
Because of the T-reversal symmetry, each eigenstate of this, say \tp n (r)), has its own 
Kramers-paired state \ip n (r)) : (ip n {r)\ = —is y \ip n (r)) . Thus, the single-body Green 
functions obey the following relation; 




This situation is reminiscent 




subjected 




(5) 



where '+/— ' sign is for the retarded/advanced Green function G R / A . 



I which corresponds to the inter- 'eigenvalue' scattering in the current context. 
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The quenched disorder average is taken at the Gaussian level with the short-ranged 
correlation; 

^^Jfj PM '" 6XP d3rA i,i v i(r)vi(r)]. (6) 

3,1 

Green functions thus averaged acquire the translational symmetry, G R ^ (r, r'; //) = 
When Fourier-transformed, they can be expanded in terms of the 16 
7-matrices with some complex-valued coefficients, 

4 42 45 

G R (k;/j,) = ^2%F n (k;/j,)+ ^ 7jFj ^ 7 m F m (A;; //) + ^ 7^(A;; /x), 

n=l .7=12,13,- m=0,5 1=15 

and G A (k;fi) = {G R (k; /i)} t . The T-symmetry requires F 1)2 ,3,4(fc) and F 12 ,i3,i4,23,34,42(A;) 
to be odd functions of k, while F 0j5 (k) and F 15,25,35,45 (&) to be even in k. Thus, only the 
latter six participate in the self-consistent Born equation, 

t R { l i) = G R '-\k- l i)-G R '-\k- l i) = Y J ^,i J dk'^j ■ G R (k'; n)-% 

3,1 

r 4 

= A 3,i / ■ {%h(k') + %F 5 (k') + ^7„ 5 F n5 (fc')} ■ % (7) 

3,1 n=l 

The bare Green function is defined as, 

3 

G R '~ l (k; fi) = (p + iS)% ~ k 3^J ~ m 7s = f j7j- 

3 1 J=0,l,-,5 

The inverse of the Green function thus determined is at most a linear combination of 
To, 75, 7i5, ■ ■ • and 745, 

4 

G R '-\fi) = F 7o + F575 + F ^' 5 - 

3=1 

Eq. (J7|) determines their coefficients, 

4 

F 



f = -|a ,o + A 5j5 + ^ A j5ij5 } • A rf 3 /c'F (A; / ) 

3=1 

-2A 0i5 / C / 3 fc'F 5 (A; , )-2^{Ao, i 5 / d 3 k'F j5 (k') 

4 . 
F 5 - f 5 = -{ A ,o + A 5 , 5 - A ^ 5 } " / d 3 k'F 5 (k') 

-2A 0i5 / rf 3 A;'F (A;')-2^{A 5 , i5 / d 3 k r F j5 (k') 
J j= i J 

4 . 
Fj5 — — I Ao,o — A5_5 — A m 5 jm 5| • / d^k'Fj^k') 

m=l ^ 



(8) 



(9) 
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-2A 0j5 / d 3 k%(k') - 2A 5>j5 I d 3 k%(k')-2j2{^ m5 , j5 I d 3 k'F m5 (k')y (10) 
That is, the right hand sides are given by these coefficients themselves, 

F„ r _ _ . 4 

F (fc) : 



<K*0 



{^ 2 -F^ + F 2 + £F 2 5 }, F B (*) 



F i6 (fc) 



"i5 



with gr(fe) being defined as 



i=i 

t 2 + F 2 -F 2 -VF' 



2 

?/(5 



m=l 



. 2 /c,i 



^ ^ Fm5^V. 



m=l 



9(k) = | - Fq + Fg + ^ F 2 5 - /c 2 | 2 + 4fc 2 | F5 - Fq| + 4| F m5 fc m } 2 - 

j=l m=l 

Note that all integrals in eqs. ( 18J9|) depend on the ultraviolet cutoff; / <i 3 /c = J| fc i <A ^ 3 &. 

We will solve these coupled integral equations, assuming that the spatial inversion 
symmetry is recovered after the quenched average. Namely, we assume 

%-G R ^ A \k^)-% = G^ A \-k^), (11) 

because only 75 anticommutes with both 71,. ..,4 and 715,. ..,45. This requires A j5, A 5j - 5 
and Fj 5 (J — 1, ■ • • , 4) to be zero, so that the equation becomes, 

1 k 2 dk 



F 
F 5 



+ G 



B 



-B 
-A., + A n 





' F ' 






, G = 2 f 








—m 


Jo 



F 2 - F 2 



k 2 



(12) 



with A s = 2vrA(A ,o + A 6 , 5 ), A a = £J =1 2vrAA j5ij 5, and B = 47rAA , 5 . In eq. (jrj), the 
momenta and energies are rescaled by the ultraviolet cut-off A, 

A ->■ 1, k -> kA' 1 , F ,5 -> ^0,5 = Fo, 5 A _1 , (/i,m) old -)> (/i,m) new = (/i, m) old A _1 . (13) 

The first part of eq. (fT2l can be 'diagonalized' by the following canonical transformation, 



" ^0 " 









ch# -shfl 
-she chfl 

where the angle 9 is chosen as 
(ch0, sh0) 





' F ' 










_F 5 _ 


5 


^2 





ch# -sh0 
-she ch9 



fi 

-m 



(A s + a/A 2 -!? 2 , fi) 



V^{A2-B2 }T{>/A 2_ B 2 + Aa} 

Under this transformation, we obtain, 



1 • 

1 



;i + r ?+ G)J- = ^i 
;i - 77_G)7- 5 = $ 



G 



2, 



fc 2 dJb 



T 2 — T 2 — h 2 



= 2 



fc 2 c/£; 



(a + *) 2 -fc 2 ' 



(14) 



(15) 



(16) 



with 77^ = a/A 2 — i? 2 ± A a and (a + ife) 2 = J 7 , 2 — The real/imaginary part of l C is 
an even/odd function of both a and b, 



ReG 



a, 

-2 - ^ 



ImG = — lot 

2 e 



;i-a) 2 + 6 2 
;i + a) 2 + 6 2 
1 - a) 2 + b 2 
:i + a) 2 + b 2 



+ 6( arctan 
1 - 



1 — ai 



arctan 



1 + a 



a ^arctan 



arctan 



1 + a 



(17) 
(18) 
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Eqs. (I16H18I) can be solved in a quite similar fashion as we did in Ref [IT] (compare 
them with eqs. (38-41) in this reference). As a result, the final self-consistent Born 
phase diagram obtained in this section is basically same as that obtained there. The 
only difference is that the phase boundaries which separate two gapped phases and 
compressible (diffusive) phase are deformed in the \i-m parameter space, according to 
the canonical transformation eq. (II 4ft: they become symmetric with respect to ipi and 
%p2i instead of /i and m (compre Fig. |5]with Figs. 9,10 in the reference). Those who are 
well-informed of Ref. [UJ therefore might as well skip the remaining part of this section 
and start from sec. 5. To make this article self-contained, we will review the arguments 
in the context of eqs. (11611181) . 

4-1. if>i — 1P2 — case 

Consider the simplest case first, jj = m = 0. Since r] + + r]_ = 2 a/A^ — B 2 ^ [§[ the 
possible solutions of eqs. ( TT6lfT8|) are clearly three-folded: 

(i) : Jo = ^5 = 0, 

(ii) : 1 + r] + G = n Js = 0, (19) 

(iii) : 1 - r]_G = n Jo = 0. 

The type- (iii) solution can realize when r/_ < — ~. In reality, however, such a parameter 
region is very limited and the solution itself is also physically unlikely. Thus, we will 
investigate only the first two types. 

The type-(i) solution is a diffusionless solution, where the massless zero-energy 
state is free from any disorders. The type-(ii) solution describes the diffusive massless 
state: the zero energy state at the critical point acquires a finite life-time due to the 
non-magnetic disorders. The gap equations have an intrinsic critical disorder strength 
T] +}C , below which only the type-(i) solution is allowed. When the disorder strength 
exceeds this critical value, r/ + > i] + c , this type-(i) solution becomes unphysical and type- 
(ii) solution becomes a unique physical solution. These two solutions are continuously 
connected at r] + = r] +)C . 

To see this situation, begin with the type-(ii) solution: 

1 + ri+ReG = n ImG = n J 5 = 0. (20) 

ImG = requires either (a) a = or (b) b = fl \a\ > 1. Now that 'a' and '6' as well as 
'F ' and 'F 5 ' are rescaled by the ultraviolet cutoff A as in eq. (113"]) . they should be much 
smaller than unit. This leads to a = 0. With this, 1 + r] + R,eG = determines '6' as a 
function of the disorder strength r) + : 

b ■ arctanfr 1 ! = 1 - — . (21) 

This equation dictates that '6' can be finite only when r] + is greater than some 
critical value. That is r) +tC = |. Above this, the type-(ii) solution becomes possible, 

§ The convergence of the gaussian integral in eq. © requires A s > B. 
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(J"0) ^5) — (ib, 0). When the disorder strength falls below this critical value, the solution 
reduces continuously to the trivial one, (To, ^5) = (0,0). In the followings, we will 
summarize the behaviours of these two solutions in the presence of finite ip\ and 1^2 ■ 

4-2. ipi ^ and ip2 — case 

Introduce small ip\ into the type-(i) solution first. Employing J-5 = 0, we have only to 
solve the first line of eq. f lT6|) in favor for a + ib = To. Since its real part, T' Q) is an odd 
function of ipi and its imaginary part J-' Q ' is even in ipi, 'a' is the first order in small 
ipi while '6' is the second order, a = 0(ipi) and b = 0(ip\). With this, eq. ( |T6l) can be 
solved iteratively in small 

This indicates that, for n + > i] +jC , the sign of the renormalized chemical potential 'a' 
and that of the bare chemical potential becomes opposite with each other. This 
is, however, physically unlikely. When rj + exceeds this critical strength, the type-(ii) 
solution instead of type-(i) becomes a physical solution. Indeed, having a finite '6' 
already at the zero-th order in ipi, the type-(ii) solution behaves as follows, 



+ O^'l fe^i), b = — * — + o(ip(, r), (23) 

2TJ+ - 1 7T T] + 

in small ipi region. For r] + > i] +)C , 'a' has the same sign as that of the bare one. 
4-3. ^i=0 and ip2 7^ case 

In the presence of finite ip2, the solutions of eqs. (TTMTBj) are two-folded, depending on 
the relative ratio between the disorder strength and the topological mass ip2- When 
the topological mass is less than a certain critical value, say ip2,c, the system is in the 
diffusive (compressible) phase, which is basically equivalent to the type-(ii) solution at 




l-^)', (24) 



r and ip2,c are defined by rj±: 

tfj 2 , c (v + , V-) = V+ + V ~ , r-Wantr] = 1 - -L (25) 
r/+r 2r/ + 

When ^2 exceeds this critical value, the system enters one of the two incompressible 
phases, which is adiabatically connected to the gapped phases in the clean limit, 

Jo = 0, F 6 = (26) 

The renormalized mass terms ip 2 * s given by 

^ 2 {1 + 2?7_ - 2?7_ ; ^ 2 arctan[7 2 ~ 1 ]} = ijj 2 . (27) 

These two solutions are connected continuously at the phase boundary, ip 2 = 
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4-4- 4>i 7^ and ip2 7^ case 



yfj>0, y/ 2 <0 

(To, To, -r;, -Ts) .retarded 
(To ,-!F*g '-Fs.fs ) - advanced 



I 



y/j<0,y/ 2 <0 .4 

(-r \r;,-r;,r;) -.ret. 
-r;, -r;, -r;) .adv.j 

I 



I 

I 

I 

-4 

/ 



/ 



/ 
/ 

r 

/ 

/ 
/ 



¥i 



Wi>0, y/ 2 >0 

/ (r ,-ro,r;,-r;) -.adv. 



¥2 



> m 



y/j<0, yr 2 >0 

(-r;,rj,r;,-r;) -.ret. 
(-^„,-r;,r;,r;') .adv. 



Figure 3. J-q = J-q + %Fq and T§ = T' b + iT'{ as a function of bare chemical potential 



H and topological mass m. changes its sign under ipi — > 
sign under ip2 — > — ip2- 



--01, while J5 changes its 



For general ipi and ip2, we have solved numerically the following coupled equations 
in favor of 'a' and '£>': 

-0i(l + rj + ReG) . ipiTj+ImG 

J~ — — : ^ / T ^TTTT + 1>- 



[1 + r] + ReG) 2 + {r] + lmG) 2 (1 + r] + ReG) 2 + {r] + lmG) 2 ' 
^2(1 + 77-ReG) . ip 2 7]_lmG 



(1 + r/_ReG) 2 + (r/_ImG) 2 ' " (1 + 7?_ReG) 2 + (r/_ImG) 2 ' 

with — J 72 = (a + i&) 2 and ReG and ImG defined in eqs. (117|ll8j) . The numerical 
solutions thus obtained are clearly four-folded, i.e. (a, b), (—a, 6), (a, —b) and (—a, —6). 
This leads to double degeneracy in and J-5 

with = F'o + zJ-"o an d ^5 = ^"5 + ^5 • Moreover, these two solutions at (^1, V^) are 
related to those at the other three points, i.e. (/0i, —1P2), ( — ^l; ^2) and (—^1, —^2), 

(J 7 3' 3 > "^5> d)|^i,^2 = 3' 3 5 — ^ 5' "F^ 7 5 ) "02 

— ( — J 7 3) 3 > ^55 5 ) I —^1,^2 = ( — -^i 3' 3 > — J 7 55 i*^ 7 5 ) I --01, --02- (28) 

The upper/lower sign is for the retarded/advanced Green function (Fig. [3]). 

Fig. H] demonstrates how Tq and J-5 behave as a function of and m, both for 
?7 + < r] +tC and for 77+ > rj +yC . Fig. [5] is the corresponding phase diagrams, where 
both J-q and T'l vanish in two incompressible phases (the topological insulator and 



disordered topological quantum critical points in three-dimensional systems 15 




-0.2 -0.1 0.0 0.1 0.2 -0.2 -0.1 0.0 0.1 0.2 -0-2 -0-1 0.0 0.1 0.2 -0.2 -0.1 0.0 0.1 0.2 

m m m m 



Figure 4. (^0,^5) numerically evaluated at A Q /A S = B/A s = 0.5 and 
774. = 0.45/0.57 (upper/lower four panels). Their real parts and imaginary parts 
decrease/increase toward the dark/light region, (upper/lower four from left to right): 
The contour plot of Tq, Tq , J 7 ^ and T' 5 ' at t] + = 0.45/0.57. The contour intervals are 
0.32/0.30 x 10-\ 0.9/1.25 x 10~ 2 , 0.24/0.24 x 10" 1 and 0.68/0.75 x 10~ 3 respectively 
Both T'q and become zero at ipi = 0. Fq and T§ vanish within both the yellow 
region (topological insulator) and the red region (ordinary insulator) . JFq in the lower 
panel remains Unite even at "01 = 0, as far as —ip2,c < 1P2 < V^.c, where ^>i,c is given 
by eq. (|2"5)). 




Figure 5. The self-consistent Born phase diagram for fixed A a /A s = B/A s = 0.5. 
(a/b/c) r/+ = 0.31/0.48/0.53, where ry+, c = 0.5. 



disordered topological quantum critical points in three-dimensional systems 



16 



an ordinary insulator). For r/ + < r) +tC , the \i-m phase diagram holds a direct phase 
transition point between these two gapped phases (Fig. Eta,b)). For r/ + > rj +;C , this 
direct phase transition point becomes a finite region of the diffusive phase (Fig. 0(c)). 
Especially at ip\ = 0, this diffusive region ranges from ip 2 = — ip2,c to ip 2 = i>2,c, where 
ip2,c is given by eq. (1251) . Note also that, due to the symmetry relation given in eq. ( 1281) . 
all the phase boundaries are symmetric with respect to the sign changes of ipi and ip2 
(see red arrows in Fig. E](c)). 



5. quantum conductivity correction in the presence of the topological-mass 
type disorder potential 

We will argue the behaviour of the quantum conductivity correction in the presence of 
both chemical-potential type disorder and topological-mass type disorder. To do this, 
consider the following series-sum of the ladder type diagrams; 

r dif (g,u;) = [r 1 -J2 dR ( k +^+) x G A {k^^)Y (29) 

k 

where k± = k ± |, jj± = /j, ± ^ and the 16 by 16 tensor | is given by [J]; 

^ = A ,o7o x 7o + A , 5 75 x 7o + A 5i0 7 x % + A 5:5 % x %. 

Tracing out its two vertices with some internal degree of freedom, say jj, we obtain the 
correlation function of the corresponding density ip^jip, 

Ml*") = ~2ri Y,Y,^^Gl(k + ,fi + )G^(k-,fi.) 

a,-,0 k,k' 

{5 k , k >6 e5 6 jC + rg )< (g ) w)^(A/ +) Ai + )^(* , _,A*-)}[7i]* y (30) 

When this density is a conserved quantity in each ensemble, like ip'joip, the function 
should exhibit the diffusive behaviour at the infrared region, q,u ~ 0. Thus, the above 
series-sum carries at least one diffusion pole. When its advanced line is T-reversed, it 
is transcribed into the series-sum of the 'maximally-crossed' diagrams (Cooperon), 

ij=0,5 

i,---,m k 

+ Y,^2 Ai >3 Al ' mAn 'P Na«i^f l5l (A; + ,/i + )[7i]5 1 a 2 Gf 252 (/c+,/i + )[7n]5 2( 5 
i,—,p k.k' 

x [7i] 77l G^ 1 (-A;_,^)[7 m ] / 3 l72 G4 & (-^,/i-)[7 P ]/3 2 /3 + •••, (31) 

which describes the quantum interference between the time-reversal paired scattering 
process from K to — K + q. Thus, any type of the diffusion pole appearing in eq. 



|| Note that the product between these tensors is defined as (A r x A a ) ■ (B r x B a ) = A r B r x B a A a , 
where A r (B r ) is a 4 by 4 matrix in the retarded line, while A a (B a ) stands for that in the advanced 
line. 
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generally results in a certain quantum interference effect in the backward scattering 
channel. 

To investigate the nature of the diffusion poles in eq. ( 129]) . employ the canonical 
transformation used in sec. 4 first, 

T di \q,co) = V* x V e A ■ T dil (q,co) ■ V e R x V e A , 

9 9 
v e R = v e A = cosn 2 7o + sinn 2 

with 9 being defined as eq. ( |T5l) . Such a transformation diagonalizes £. Simultaneously, 
it simplifies the expression of G R ^ (k±, /j±) in terms of J-'o and J- 5 . Namely, the series- 
sum thus transformed reads 

T^M = -Y^G R (K,» + ) xG^k-^-)]^ (32) 

k 

with 

€d = — 2 70 X 70 2 2 75 X 75, 

3 

r>Rn \— •'0,+ * v - ^ 7 - ^ r 5-+ 
y {k+, ft+) = 2 _ , 2 _ t2 7o - 2^ fc +-i7j - ^5 ^75, 



T 2 — h 2 — T 2 ' u /-^ '"~ r,J ' J 7^ 2 — _ T 2 
- ^+ ^5,+ j =1 - r 0,+ •'5,4- 

X"* 3 x> 



0A{k - ^ ~ (^_) 2 - k[- (^_) 2 ^° g fc "^ (^_) 2 - (^_) 
2i; ± = VA 2 - 5 2 ± 27t(A 0i0 - A 5)5 ). 



The subscript on stands for its w-dependence; T^± = ± %,m), where the /i- 
and m-dependence of J-j are determined by the preceding gap equation, eqs. (116til8j) . 
Taking q to be zero, we get a relatively simpler expression for eq. fl32l) . 

f dif (0, oo) = /r 1 ^ + / 2 - x r 2 + /g-^a + f^T i: (33) 

/• 2 r 2 , r 2 -f_2 2,2 

/i — a x — oa 04 -+- oa 23 , j 2 — a 1 — a 04 -t- a 23 , 
/ 3 = 9a 2 - 5^4 + <5 a L h = 9a 2 - + a 23> 
with 5dij ee <2j — a,-, ajj ee + a,- and 

_ ^+ _ \^ A+^o- r o 4 x 

a °-<-^ 2 ^(f-j- 2 + + j- 2 + )(P-(j- v) 2 + (J- 5 v) 2 )' 1 } 

k 2 

^ = " V (* 2 - ^ + ^K*^ (^O,-) 2 + (^5,-) 2 ) ' ^ 

Q2 = ^ (* 2 - ^0% + - (^O,-) 2 + (^5,-) 2 ) ' (36) 

J~5 J~* 

^ ^ ^ (* 2 - ^0 2 + + ^, + )(* 2 " (^0,-) 2 + (^5,-) 2 ) ' 

v- v-^ ^ r 5,+-7 r 5 



a 4 = 



_ V ,+ 5 '~ f3£ 

^ fP _ r 2 -i- T 2 \(h 2 -(T* Vm(T* V»v v 
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Four Tj in eq. (13"3"j) turns out to be all regular functions (tensors) of uj [llj. The only 
possibility is, therefore, the diffusion pole appears in one of the four coefficients fj. In 
fact, employing the same technique invented in Ref. [TT], we can easily show that eq. (1TB1) 
requires to be zero at uj = 0. The other three coefficients are generally truncated by 
some finite infrared cutoff. Each cutoff stands for (the inverse of) the relaxation time 
of a certain internal degree of freedom other than charge density. 

When the topological-mass-type disorder is absent [TTJ , f 3 also exhibits the diffusive 
behaviour at the massless point, m = 0. This is because, in the presence of only the 
7 -type disorder potential, 745 always commutes with a hamiltonian at the TQCP, so 
that not only the charge density but the corresponding parity density^ also follow a 
diffusion equation. Indeed, when substituted into eq. fl30|) . T 3 contributes to the 
correlation function of this parity density, while / 4 ~ 1 r 4 does to that of the charge 
density. When their advanced lines are time- reversed as in eq. ( )3T|) . both T 3 and 
result in the same magnitude of the positive weak-localization (AWL) correction to the 
charge conductivity. Thus, the additional U(l) symmetry recovery at the TQCP induces 
the 'doubling phenomena' of the quantum conductivity correction around this massless 
point. 

In the presence of the topological-mass-type disorder, however, the relaxation time 
of this parity density, given below, always remains finite in the whole fi-m parameter 
space, 

<o - {^(/a - /4)}^ o = {^(-«o«4 + ^3)}^. (39) 

Namely, though both a 2 and 03 in the r.h.s. vanish at ip 2 = — /i sinh # + m cosh # = 0, 04 
never vanishes even at ip 2 = 0, due to the first term in eq. (138|) . — -rzrn ■ Indeed, once the 
75-type impurity potential is introduced, the parity density is generally a non-conserved 
quantity. As a result, unlike in Ref. only the charge diffusion pole results in the 
positive weak-localization correction, while that of the parity diffusion mode remains 
always 'gapped' in the entire /i-m parameter space. 

6. summary 

From the bulk-edge correspondence, the topological quantum critical point (TQCP) 
which intervenes the three-dimensional topological insulator and an ordinary insulator 
is expected to be stable against any non-magnetic disorders, provided that each 
surface state supported in the topological insulator phase is stable against the same 
perturbations. To understand the stability of this 3-d TQCP, we first employ the 'Berry 
phase' argument and show that any backward scattering process which conserves the 
parity density degree of freedom, ip^^ip, is always set off by its T-reversal counter- 
process. This observation upholds more directly our previous surmise in Ref [UJ, where 
two of the authors conjectured that, when the system transits from the topological 

745 is spatially-inversion odd, while time-reversal even, so that we call the corresponding density, 
^ 745^1 as the parity density. 
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insulator side to the ordinary insulator side, the parity density always becomes a 
conserved quantity once. However, it is still an open issue how the parity density 
becomes a conserved quantity at these transition points (or regions) in the presence of 
generic non-magnetic disorders. Namely, the 'Berry phase' argument only confirms that, 
as far as the parity density is a conserved quantity, the 3-d TQCP remains delocalized 
(or critical) even in the presence of sufficiently strong generic non-magnetic disorders. 

In fact, as for those backward scattering processes which do not conserve the parity 
density, the above 'Berry phase' argument does not work at all. To understand these 
generic T-symmetric situations, we further derive the self-consistent Born phase diagram 
in the presence of general non-magnetic disorders. We found that the derived scB phase 
diagram is basically same as that in the case with only the chemical-potential-type 
disorder. Namely, as in Ref. [11] . there exists a certain critical disorder strength, below 
which the two gapped phases (the topological insulator and an ordinary insulator) are 
always separated by the direct quantum phase transition point in the \i-m parameter 
space. When the disorder strength exceeds this critical value, this direct transition point 
becomes a finite region of the diffusive metallic phase. 

As for the quantum conductivity correction, the situation changes drastically. On 
increasing the topological-mass type disorder potential, we found that the diffusive 
parity density mode observed at the 3-d TQCP always acquires a finite relaxation time 
in the entire fi-m parameter space. This indicates that the 'doubling phenomena' of 
the quantum conductivity correction previously observed at the TQCP [11] becomes 
suppressed in the presence of 75-type disorder potential. 
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